Kaon 5-parameter from improved staggered fermions in Nf = 2 + 1 QCD 
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We present a calculation of the kaon B-parameter, Bk, using lattice QCD. We use improved 
staggered valence and sea fermions, the latter generated by the MILC collaboration with Nj =2 + 1 
light flavors. To control discretization errors, we use four different lattice spacings ranging down 
to a ~ 0.045 fm. The chiral and continuum extrapolations are done using SU(2) staggered chiral 
perturbation theory. Our final result is I3k ~ 0.727 ± 0.004(stat) ± 0.038(sys), where the dominant 
systematic error is from our use of truncated (one-loop) matching factors. 



CP violation was first observed in the kaon system in 
1964 1] , but only in the last few years has it been possible 
to use this classic result to learn about the parameters of 
the Standard Model (SM). CP violation in the SM is due 
to the phase in the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix, and leads to CP violation in kaon mixing — so- 
called indirect CP violation — through Feynman diagrams 
involving virtual charm and top quarks. Integrating out 
heavy quarks (c, b and t) and the W and Z-bosons, one 
finds that a prediction for indirect CP violation requires 
the calculation of the matrix element of a local four- 
fcrmion operator between a Kq and a Kq. To calcu- 
late this matrix element (which is parametrized by the 
kaon i3-parameter, Bx) requires control over the non- 
perturbative physics of the strong interactions. The only 
known quantitative method to calculate such matrix ele- 
ments from first principles is lattice Quantum Chromody- 
namics (LQCD), and it is only very recently that lattice 
calculations have begun to control all sources of error. 

In this note we present a calculation of Bk using im- 
proved staggered fermions. Fully controlled results us- 
ing other types of lattice fermion have been obtained 
previously, using valence domain-wall quarks with stag- 
gered sea quarks [H, using domain- wall valence and sea 
fermions 0,3. and, very recently, using Wilson valence 
and sea fermions Q. We also presented a result using a 
partial data set in Ref. @, which we refer to as SW-1 in 
the following. What sets our work apart from that using 
valence domain-wall fermions is that we are able to use 
smaller values of the lattice spacing, a, thus potentially 
providing better control over the continuum limit. We 
also have very small statistical errors, such that the final 
statistical error in Bk is ^ 0.5%. In addition, our use of 
a different fermion discretization provides a highly non- 
trivial cross-check of the other results, analogous to the 
use of a different experimental technique. 



The required matrix element is parametrized by 

~ 8/iA4/3 

where R is a specific regularization scheme chosen to de- 
fine the operator Oas=2 = HA^lA^ ~ lz)d][s^v{l - 
75)(i], and is the corresponding renormalization scale. 
We use lattice regularization, and then convert to the 
continuum MS scheme (using naive dimensional regular- 
ization for the 75) using one-loop matching factors from 
Ref. At the end we convert to the renormalization- 
scale invariant quantity Bk- 

We use improved staggered fermions for both valence 
and sea quarks. The advantages of staggered fermions are 
that they are computationally inexpensive and that they 
retain a remnant of chiral symmetry. The latter prop- 
erty implies that the matrix element in the numerator of 
Eq. ^ vanishes when Mk because of the "left-left" 
chiral structure of Oas=2- Without chiral symmetry, 
Oas=2 mixes with operators with "left-right" structure, 
whose matrix elements do not vanish when Mk 0, and 
are thus enhanced. For both staggered and domain-wall 
fermions, such mixing is not allowed. Mixing with chi- 
rally enhanced operators is allowed for Wilson fermions, 
but appears now to be controllable [s]. 

For the sea quarks, we use the MILC collaboration's 
publicly available ensembles generated using 2 + 1 flavors 
of asqtad staggered fermions Q. Here "2 + 1" indicates 
degenerate up and down quarks and a heavier, nearly 
physical, strange quark. Each staggered lattice flavor de- 
scribes four continuum fermions, usually called "tastes" . 
This unwanted degeneracy is removed by the fourth-root 
prescription, and we assume this leads to no problems 
with the continuum limit. 

For the valence sector we choose HYP-smeared stag- 
gered fermions. We prefer these to asqtad fermions be- 
cause they are more continuum-like, e.g. the breaking of 
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TABLE I. MILC ensembles used to calculate Bk iSj. a is the 
nominal lattice spacing, me (ms) the light (strange) sea-quark 
mass, "ens" the number of gauge configurations and "meas" 
the number of measurements per configuration. "Status" in- 
dicates changes since SW-1: "old" is unchanged, "new" is a 
new ensemble, and update indicates more measurements. 



a (fm) 


ame/ arris 


size 


ID 


ens X meas 


status 


0.12 


0.03/0.05 


20^ X 64 


CI 


564 X 9 


update 


0.12 


0.02/0.05 


20^ X 64 


C2 


486 X 9 


update 


0.12 


0.01/0.05 


20^ X 64 


C3 


671 X 9 


old 


0.12 


0.01/0.05 


28^ X 64 


C3-2 


275 X 8 


old 


0.12 


0.007/0.05 


20^ X 64 


C4 


651 X 10 


old 


0.12 


0.005/0.05 


24^ X 64 


C5 


509 X 9 


update 


0.09 


0.0062/0.031 


28'' X 96 


Fl 


995 X 9 


update 


0.06 


0.0036/0.018 


48^ X 144 


SI 


744 X 2 


old 


0.045 


0.0028/0.014 


64^ X 192 


Ul 


705 X 1 


new 



taste symmetry is smaller by a factor of ~ 3 Nev- 
ertheless, taste-breaking induced by discretization errors 
leads to significant complications in the analysis. This 
enters both through mass splittings between kaons of dif- 
ferent tastes and by inducing additional operator mixing. 
Both these effects, as well as those due to the use of dif- 
ferent types of valence and sea quarks, and to the use of 
the fourth-root prescription, can be incorporated into the 
chiral effective theory describing staggered fermions — 
staggered chiral perturbation theory (SChPT) (gI-IToI- 12[. 
Specifically, we use SU(2) SChPT — in which only the up 
and down quarks are treated as light. After fitting to 
the forms predicted by SChPT, one can then remove "by 
hand" the taste-breaking discretization errors. 

We use the MILC asqtad lattices listed in Table U for 
the present work. The most important changes since our 
previous work, SW-1, are the addition of a fourth, finer, 
lattice spacing (the "ultrafine" ensemble Ul) and the 9 
or 10-fold increase in the number of measurements on 
several of the "coarse" ensembles (a « 0.12 fm) and also 
on the "fine" (a « 0.09 fm) ensemble Fl. 

The sea-quark masses in these ensembles are not phys- 
ical. The strange quark is somewhat too heavy, re- 
quiring a small correction. The light quarks are too 
heavy, requiring an extrapolation to the physical mass. 
The lightest sea-quark pion on the coarse ensembles has 
m™'"(sea) « 280 MeV. This will turn out to be light 
enough for a controlled extrapolation, because the de- 
pendence on sea-quark masses turns out to be mild. 

On each lattice, we create and destroy kaons using two 
wall-sources, which are separated by Euclidean time At. 
The sources have the property of creating only kaons with 
the desired Goldstone taste, ^5, and with vanishing spa- 
tial momenta. The discretized version of the operator 
Oas=2 is placed between the two sources, and summed 
over space. At is chosen large enough that there is a 
plateau region with minimal contamination from excited 
states, and small enough that effects from states propa- 



gating "around the world" in the time direction can be 
ignored. We obtain Bk by fitting to a constant over the 
plateau region. 

We use 10 different valence quark masses, amv — 
ar7i"°™(n/10) for n = 1,2,..., 10, with the nominal 
strange quark masses being 0.05, 0.03, 0.018 and 0.014 
for coarse, fine, superfine and ultrafine lattices, respec- 
tively. (For equal bare masses, HYP-smeared quarks 
have smaller physical masses than asqtad quarks, because 
of differing renormalization factors.) We use the lightest 
four valence masses for the valence d quark (nix), and 
the heaviest three valence masses for the valence s quark 
{rriy). This maintains the relations rrix ^ my ^ m'^^'^^^ 
as required for the applicability of SU(2) ChPT. 

The four values of aurix allow us to extrapolate to the 
physical down-quark mass. This extrapolation is much 
shorter than that for the sea-quarks, both since our va- 
lence quarks are lighter (the lightest xx pion has a mass 
« 200 MeV) , and because we are extrapolating to mF^^^ 
rather than {m^}^^^ +mF^^^) /2 (so that M^x must be ex- 
trapolated to 158 MeV, the mass of an unphysical dd 
meson). We fit the dependence on X = (for pseu- 
doscalar taste) to the next-to-leading order (NLO) form 
predicted by SU(2) SChPT. In SW-1 we used uncorre- 
lated fits and did not include finite-volume corrections. 
Here we correct both shortcomings. To obtain satisfac- 
tory correlated fits (with x^/d.o.f. < 1) we need to in- 
clude higher-order terms with coefficients constrained by 
Bayesian priors. Specifically, we fit Bk for fixed ruy to 

CiFo(X)+C2A:+C3A:2+c4A:2 In^ X+c^X^ \nX+CfiX^ , (2) 

where Fq{X) contains the leading order constant term as 
well as the chiral logarithms. The latter include taste- 
breaking effects and finite- volume dependence (see SW-1 
and Ref. [l3| for the explicit form). The terms multi- 
plied by C3_5 are the generic NNLO forms in the contin- 
uum. Since these are not known analytically, we include 
them with coefficients whose magnitude is constrained 
by Bayesian priors to be of the size expected by naive 
dimensional analysis. We also include a single analytic 
NNNLO term (with coefficient Ce). 

Examples of the resulting fits are shown in Fig.[TJ With 
the fit parameters in hand, we can extrapolate rrix and 
rrf^'^ to mF^^^, rrf°^ to ml^^^, the volume to infinity, and 
remove taste-breaking discretization errors. We note that 
finite volume shifts are small for our pion masses, even 
though our lightest pions have to^L « 3 (with L the box 
size) . We have checked this directly by comparing results 
on ensembles C3 and C3-2 [l^. 

We estimate the systematic error of our X-fits by dou- 
bling the allowed widths of the Bayesian priors, and by 
dropping the NNNLO term. The former gives the larger 
effect, and we take the largest size of this shift (0.33%, 
on the SI ensemble) as the error estimate. 

The three values of amy allow us to extrapolate to the 
physical strange quark mass. We find that a linear fit 
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FIG. 1. _Bk(NDR, 1/a GeV) vs. X on the Ul and Fl 
ensembles. The [red] diamond ([blue] cross) shows the result 
on the Ul (Fl) ensemble, after extrapolation and removal of 
taste-breaking artifacts as described in the text. 



to my represents the data very well, and use this for our 
central values. We use a quadratic "Y-fit" to estimate a 
fitting systematic. 

At this stage the values of Bk on different ensembles 
differ primarily because of taste-conserving discretization 
and matching errors. To remove the main part of these 
errors we use ensembles C3, Fl, SI and Ul, which have 
very similar sea-quark masses. In Fig. [21 we show the de- 
pendence on a^, and present the results of several meth- 
ods of extrapolation to a = 0. We note that the simple 
linear dependence observed in SW-1 (for the largest three 
lattice spacings) has been resolved by improved errors, 
and by the addition of the Ul point, into a less smooth 
dependence. On theoretical grounds 
dependence is 



12|, the expected 



di + d2{aKf + d^iaAfas + d^al + dz{aKf + . . . , (3) 

where — a^^(I/a). We fit to this 5-parameter form 
applying Bayesian constraints on 1^2-5 — the expected val- 
ues are taken to be 0, while the standard deviations are 
set to 2, having chosen A = 300 MeV. The fit is shown 
by the [blue] dotted curve, and gives the extrapolated 
value shown by the diamond. The fit, however, is very 
poor, with x^/d.o.f. — 4.5. This problem is not resolved 
by adding terms of one higher order. Thus we drop the 
coarse lattice from the fits, and then find good fits to 
a constant (solid [red] line and cross), a linear depen- 
dence on a? (not shown) and to the constrained form 
([3]) (dashed [brown] line and square). We take the con- 
stant fit for the central value, and the difference between 
it and the constrained fit as the systematic error in the 
continuum extrapolation. For more discussion of fits see 
Ref. 
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After the preceding analysis, Bk can still have a resid- 
ual dependence on the sea quark masses. In SU(2) 



FIG. 2. Continuum extrapolation of Ba'(NDR,2 GeV). The 
fits are described in the text. 



SChPT, the dependence on mi is linear at NLO. We have 
investigated the mi dependence in detail on the coarse 
lattices, with results shown in Fig. [3] We plot versus 
Lp, the squared mass of the sea-quark pion, and find a 
linear behavior with a small slope ~ —1/(2.9 GeV)^. In 
SW-1, with errors 3 times larger, we could not uncover 
this dependence. Using this slope, we find that Bk is 
increased by 1.5% when changing Lp from its value on 
ensemble C3 to its physical value. Since this is a small 
effect, and since we only have results for the slope on the 
coarse ensembles, we do not adjust the central value of 
Bk, but instead quote the 1.5% as a systematic error. 

We also need to correct for the mismatch between ams 
and am^^^^. Here we follow SW-1, and assume the SU(3) 
ChPT form of the mg dependence. Then we find that, on 
C3, the correction from using a sea-quark mass which is 
40% too large is 1.3%. This approach is more conserva- 
tive than that used in SW-1. Again, we do not adjust the 
central value, but quote the shift as a systematic error. 

We collect all our errors in Table HIl The statistical er- 
ror has been reduced by a factor of 3 compared to SW-1, 
and now is smaller than many other sources of error. The 
dominant error comes from our use of one-loop matching. 
We estimate this as SBk/Bk = a^, with evaluated 
at the scale of our finest lattice. It is thus reduced from 
the 5.5% estimate of SW-1 by the addition of the Ul 
ensemble. The discretization error is nearly unchanged 
from SW-1, despite the addition of the ultrafine lattices, 
because of the increased uncertainty in the continuum 
extrapolation. The only errors not discussed above are 
the "ri" and "/tt" errors, which are estimated essentially 
as discussed in SW-1. 

Adding systematic errors in quadrature, we find 

Bk = 0.727 ± 0.004(stat) ± 0.038(sys) . 

The central value is almost unchanged from SW-1, but 
the significant improvements we have made have both re- 
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FIG. 3. i?K(NDR, 2 GeV) vs. Lp on the coarse ensembles. 
The solid line shows a linear fit, with the extrapolation to 
physical pion mass given by the diamond. 
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FIG. 4. Comparison of our result for Bk ~ -Ba'(RGI) with 
other results obtained with Nf =2 + 1 flavors. 



TABLE II. Error budget for Bk using SU(2) SChPT fitting. 



statistics 

matching factor 

discretization 

X-fits 

Y-fits 

anil extrap 

anis extrap 

finite volume 

ri 



(%) 



memo 



0.6 see text 

4.4 Asi?' (Ul) 

1.9 difl. of constant and constrained fits 

0.33 varying Bayesian priors (SI) 

0.07 difl. of hnear and quadratic (C3) 

1.5 diff. of (C3) and linear extrap 
1.3 diff. of (C3) and linear extrap 
0.5 difl. of F = oo fit and FV fit 
0.14 ri error propagation (C3) 
HI 132 MeV vs. 124.4 MeV 
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ducdid the error and solidihcd our error estimates. Our 



result is consistent with other TV/ = 2 + 1 results 2, Jk, , 
as shown in Fig. |4l The largest difference is that our re- 
sult lies 1.4(7 below that of Ref ^ (the most accurate 
result). This consistency with results obtained using dif- 
ferent fermion discretizations is our most significant con- 
clusion. It is important, however, to further reduce errors 



in a [l lattice calculations to check that this consistency 



[5] 
[6] 
M 



holds up. Work is in progress to reduce our dominant sys- 
tematic using two-loop matching and non-perturbative 
renormalization. 

Lattice results for Bk now allow the venerable experi- 



mental result for en to be used to constrain the parame- 
ters of the SM. Indeed, we have now reached the situation 
that the accuracy of Bk calculations is such that errors 
from other sources dominate — in particular those from 
Vcb and the Wilson coefficient 7jcc (see, e.g., Ref. [l5|). 
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NRF grant funded by the Korean government (MEST). 



[9: 

JIO 



[11 

[12; 

[13 
[14 

[is: 



J. Christenson, J. Cronin, V . Fitch, and R. Turlay, Phys 



Rev. Lett., 13, 138 (1964) 



Rev. 



C. Aubin, J. Laiho, and R. V. de Water, Phys 
D81, 014507 (2010), arXiv : 0905. 3947 fhep-latT 

D. J. A ntonio et al. (RBC), |Ph ys. Rev. Lett., 100, 032001 
arXiv:hep- ph/0702042' 

Toki et al., |Phys. Rev., D84, 014503 (2011) 



(2008) 



Y, 



|arXiv:1012.4178 [hep -lat] 

Durr ei al. , |Phys.Lett. 



larXiv: 1106.3230 [hep-lat]' 
T. Bae et al, Phys. Rev. 
arXiv:1008.5179 [hep-lat] 
J. Kim. W. L ee, and S. R. 



B705, 477 (2011) 



D82, 114509 (2010), 



Sharpe, |Phys. Rev., D83 



094503 (2011)1 , a rXiv:1102. 1774 [hep-lat[ 
A. Bazavov et al ., Rev. JVlod. Phys., 82, 1349 (2010), 
arXiv:0903.3598 [h ep-lat] 
T. Bae et al, |Phys. Rev. 



arXiv:0801.3000 [hep-lat] 
W. Le e and S. R. Sharpe, 



D77, 094508 (2008) 



|Phys. Rev., D66, 114501 



(2002)1 , arXiv:hep-lat/0208018-" 

C. Au bin and C. Berna rd, Phys. Rev. D, 68, 034014 
hep-lat/0304014"] 
V. de Water and S. Sharpe, Phys. Rev. D, 73, 



(2003), 
R. S. 



014003 (2006), ar Xiv:hep-lat/0507012 
J. Kim et al., |Phys. Rev., D83, 



117501 (2011) 



,arXiv:1101.2685 [hep-latj. 
'W. Lee et al, PoS, "{ Lattice 2011), 316 (2011), 

'arXiv:1110.2576 [hep-lat]'! 

J. Brod and M. Gorbahn, ,arXiv:1108.2036 [hep-ph]| 



